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Introduction 


This unit is about vibrations and it continues from where Unit 8 left off. Its aim is 
to extend your ability to model systems in which vibrations play an important 
role. It will use some of the mathematical techniques you have learnt since you 
last studied vibrations in Unit 8, such as systems of differential equations, matrices 
and eigenvalues. | 


Study guide 


Each section of this unit builds on those that have gone before. You should 
therefore read the unit in the order in which it is written. The first four sections all 
contain new material, while the fifth section pulls this material together by means 
of problems which are intended to revise all the main points. You should do as 
many of these as you can (or have time for) without first looking at the solutions, 
but you should then compare your answers with those given at the end of the unit. 
Finally look at the solutions to any of the problems which you did not attempt, 
because you will probably learn something from them, too. 


There is no tape associated with this unit but there is a television programme. You 
should, if possible, have read at least the first two sections before you watch it— 
but don’t miss it even if you haven’t got that far. The programme contains 
demonstration experiments which are designed to make it easier to follow the 
theory in the unit. 


1 Modelling vibrating systems 


1.1 The story so far 


In Units 7 and 8 you studied the behaviour, both ‘free’ and ‘forced’, of some 
systems which were modelled by perfect springs, perfect dashpots and particles. 
The assembly in Figure 1 is typical of the models discussed in those units. What 
you did, in terms of Figure 1, was to write down a differential equation to 
represent the system with or without an applied force and then to obtain a 
solution in the form of an equation for x in terms of time and in terms of the 
parameters of the model (assumed constant): the stiffness k, the mass m and the 
damping coefficient r. You met a number of different models of this sort but in all 
cases the equation of the motion derived from the model took the form of a linear 
second-order differential equation with constant coefficients. You may well find it 
useful to revise some of this work by looking back at Units 7 and 8. Exercise 1 1s 
designed to help you recall some of the material in these units and to prepare you 
for later sections of this unit. If you are in any doubt about any part of this 
material you should consult Units 7 and 8. 


Exercise 1 (Revision) 


- Suppose that in Figure 1, m = 0.1;k =2 x 10° and r = 0. Write down the equation of 


motion for free vibrations and work out the angular frequency of the vibrations. 


[Solution on p. 31 | 


In Unit 7 the general solution of the equation of motion for free vibrations of the 
model in Exercise 1 was shown to be x = Xp + Csinwt + Dcoswt, where Xo, C, D 
and w are constants. (w is the ‘natural’ angular frequency which you worked out 
in Exercise 1. The constants C and D are arbitrary.) This is not the only form for 
the general solution. By choosing A and @ such that C = Acos¢@ and D= A sin 
we obtain 


Csinwt + Dcoswt = Acos ¢sin wt + Asin d cos wt 
= Asin(wt + >), 
so the general solution can be written in the alternative form 


X =X. + Asin(wt + ¢). (1) 


The behaviour of a vibrating 
system is forced if it is 
brought about by an external 
forcing agency (see Unit 8, 
Section 2); otherwise the 
behaviour is said to be free. 


Figure | 


Unless otherwise stated, all 
quantities in this unit are 
expressed in terms of SI units. 
Thus the mass of the particle 
is 0.1 kilograms and the 
stiffness of the spring is 

2 x 10° newtons per metre. 


It is usual to choose A to be 
positive and —xn<@<7,1n 
which case 

AJC +i 
and 


C 
—| ifD20 
arccos (J 1 


p = 
— arccos (J if D < 0. 
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This is a more convenient formulation for our present purpose and we shall use it 
throughout this unit. 


We can consider Solution (1) to be made up of two parts. The sinusoidal term 
Asin(at + ¢) represents the vibration. The constant Xo gives the static position 
which the particle occupies when it is not vibrating (mathematically the static 
solution corresponds to setting the arbitrary constant A to zero). The value of xo 
depends on our choice of fixed origin from which x is measured. If x is measured 
from the anchored end of the spring (as in Figure 1), xo is equal to the spring’s 
natural length. In this unit however we shall usually choose the origin to be at the 
static position so that x, = 0, thereby obtaining a simpler expression for the 
general solution. I shall have more to say about this in Section 2. 


1.2 Lumped-parameter models 


The perfect springs, perfect dashpots and particles which make up the models in 
Units 7 and 8 are idealized components. Each component represents a particular ~ 
physical property. For example, perfect springs are assumed only to have resilience 
(i.e. springiness): they deform under load in accordance with Hooke’s law and 
return to their original length when the load is removed. They are assumed to 
have no mass and no damping properties. Perfect dashpots provide linear damping 
but no mass or resilience and particles have mass but no resilience or damping 
properties. These are not, of course, completely accurate descriptions of real 
components: a real spring must have some mass and its material will exhibit at 
least some damping properties; no real body can behave exactly like a particle; it 
is bound to deform to some extent under the action of the forces to which it is 
subject—and so on. Nevertheless, models consisting of idealized components can 
be used to describe a wide variety of real systems quite adequately. 


upper torso 


mass of em 
ii et shoulder { 
system thorax- 
abdomen 
system 
springs and stiff elasticity a 
shock absorbers of spinal 
column ; 
hips 
Po mass of wheels Ered aniied 
| to sitting 
subject 
stiffness of tyres 
legs 
road 
surface annacat 
x part of tyres 
feet 
Zi force applied to 
standing subject 
Figure 2 Figure 3 


As an example of how this can be done consider the vertical vibration of a motor 
car as it moves over an uneven road surface. Figure 2 shows a deliberately 
simplified model for estimating the magnitude of these vibrations for a particular 
road profile and car speed. The model consists of idealized components each of 
which represents a physical property of the car. To keep the model as simple as 
possible we confine our attention to those properties of the car that can be 
expected to have a significant influence on its vertical motion. We can, for instance, 
expect the car’s suspension to play an important role in contributing both 
resilience and damping to the car’s vertical motion and this is represented in the 
model by a perfect spring and a perfect dashpot. Similarly, the resilience of the 
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tyres is important and is represented in the model by a second perfect spring. The 
particles in the model are used to represent the masses of the car body, the wheels, 
and the tyres. Since the model is only intended to yield information about the 
vertical motion of the car the particles are considered to move along one vertical 
straight line. 


It is important to realize that in constructing such models we do not attempt to 
describe the detailed structure of the device we are modelling. Rather we try to 
identify its important properties (mass, resilience, damping) and ‘lump’ them into 
appropriate idealized components. This kind of model is therefore known as a 
lumped-parameter model. 


The applications of lumped-parameter models are many and varied. I shall 
mention just one more example that was built up in the course of investigations 
into the effect of vibration on the human body. It was found that different 
combinations of amplitude and frequency have different effects, ranging from 
painful to imperceptible. The model in Figure 3 was found useful in predicting the 
reaction of a standing or sitting human being to vibrations at low frequencies. 
Once again it is valid only for up and down vibrations so that particles in the 
model are free to move only along vertical straight lines. 


The accuracy of the models in Figures 2 and 3 could well have been increased by 
incorporating more elements. Unfortunately this would also increase the 
complexity of the resulting equations of motion. In general, when designing 
lumped-parameter models a balance must be struck between the accuracy of the 
model and the complexity of the mathematics. 


1.3 Degrees of freedom 


Consider again the motion of the particle in Figure 1. This can be represented by 
the variation with time of one variable, x. One variable, therefore, is enough to 
specify the position of the particle at any instant: such a model is said to have one 
degree of freedom. Now look at the particle in Figure 4. Suppose that it is free to 
move in the plane of the diagram: that is to say it can move up and down and 
sideways but it cannot move into or out of the paper. One way, therefore, of 
specifying the position of the particle at any instant is to state the values of the 
displacements x and y as shown in Figure 4(a). Another possibility would be to 
state the length r of the line OP and the angle 6 as in Figure 4(b). Either way, two 
numbers (or ‘coordinates’) are required to specify the position of the particle and 
so we say the particle has two degrees of freedom. The same terminology is used 
when there is more than one particle. For example, the two spring-connected 
particles in Figure 5 are constrained to move along the straight line which 
connects their centres. So the only things that can change are the distances of the 
two particles from some fixed origin. The configuration of the system can therefore 
be specified by the two distances x, and x, shown in Figure 5. 


Exercise 2 


Make a sketch to show a different way of specifying the configuration of the two-particle 
system in Figure 5. How many degrees of freedom has this system? 


[Solution on p. 31] . 


What I have said so far can be summarized by defining the number of degrees of 
freedom of a system to be the smallest number of coordinates which are required 
to specify the configuration of the system at any instant. The importance of the 
number of degrees of freedom is that it is equal to the number of equations of 
motion. 


If a lumped-parameter model is such that each particle is free to move in just one 
direction, the number of degrees of freedom is equal to the number of particles. 
This is because the configuration of the model can be specified by giving the 
displacement of each particle along its direction of motion. For example, the 
lumped-parameter model in Figure 2 (modelling the vertical vibrations of a car) 
has three degrees of freedom because it is made up of three particles that are free 
to move vertically. 


(b) 
Figure 4 


You have already met the 
system in Figure 5 in Section 
4 of Unit 7. 


fixed line 
ee Xx] 
4 


N 


Figure 5 
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Exercise 3 
How many degrees of freedom has the model in Figure 3? 


[Solution on p. 31 | 


Exercise 4 


State the number of degrees of freedom for each of the lumped-parameter models shown in 
Figure 6 below. (The particles are constrained to move only along the lines connecting the 
particles.) 


rigid connection 


= pee 
\ (g) 
(d) tee op 


Figure 6 


[Solution on p. 31] 


Exercise 5 | 
A particle is free to move in space. How many degrees of freedom does it have? 


(Solution on p. 31 | 


The rest of this unit is an introduction to the mechanics of lumped-parameter 
models with more than one degree of freedom. We shall confine our attention to 
models of systems in which damping does not play an important role. We begin, 
in the next section, by considering a model with two degrees of freedom. 


Summary of Section 1 


This unit models real systems by combinations of idealized components each of 
which has only one of the relevant properties, i.e. mass, stiffness or damping. Such 
combinations are called lumped-parameter models. This unit deals with the analysis 
of undamped lumped-parameter models. 


The least number of coordinates that is required to define the configuration of 
such a model is called the number of degrees of freedom. 


2 Free undamped vibrations with two 


degrees of freedom particle A _ particle B 
\ 
kp 
The vibration of lumped-parameter models with more than one degree of freedom 
gives rise to phenomena which are not observed in a one-degree-of-freedom model. . 
’ This section tackles the analysis of the simplest model that will unequivocally Figure 1 


show up these phenomena. This is a model with two degrees of freedom (it is 
shown in Figure 1) and it has all the important features which are found in more 
complicated undamped lumped-parameter models. 


2.1 Setting up the equations of motion 


Before we can set up the equations of motion for the model in Figure 1 we must 
decide on the coordinates we are going to use to measure the positions of the two \ 
particles. To this end we first consider the simpler model in Figure 2 which only Figure 2 
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has one degree of freedom. If the particle in Figure 2 is displaced from its 
equilibrium position and released, it will move with simple harmonic motion about 
the equilibrium position. As I remarked in Section 1, we can obtain a simple 
expression for the motion by measuring the displacement x from the equilibrium 
position. The expression we obtain is then 


x = Asin(awt +.¢) 


where w is the natural undamped angular frequency and A and @¢ are constants 
with A positive and —z < @ < 72. (A is the amplitude — it gives the maximum 
displacement of the particle from the equilibrium position.) We could, of course, 
measure the displacement of the particle from any fixed datum point, but we 
would then obtain the more complicated expression 


X =X) + Asin(wt + ¢) 


where Xq 1s the distance of the datum point from the equilibrium position. If we 
want the simplest expression for the motion of the particle, we must measure the 
displacement from the equilibrium position. This is also true for models with more 
than one degree of freedom — indeed the simplification is even greater for such 
models, and that is why in this unit I shall always measure the displacement of a 
particle from its equilibrium position. In order to be able to interpret such 
measurements you will need to be able to find the equilibrium configuration of 
lumped-parameter models. Here are two exercises to give you practice at this. 


Exercise 1 


Each of the two springs of the lumped-parameter model in Figure 3 has a natural length of 
0.6m and a stiffness of 75 Nm’. How far are the particles from the fixed datum line when 
both particles are in equilibrium? (Remember that particles are considered to have 
negligible size.) ’ 


[Solution on p. 31 | 


Exercise 2 


The model in Figure 3 is suspended vertically from a fixed support, so that the 0.2kg 
particle is below the 0.1 kg particle and both are free to move in a vertical straight line. 
How far are the particles below the support when both particles are in equilibrium? 


[Solution on p. 31 | 


We shall now return to the model in Figure 1. This has been redrawn in Figure 4. 
Two configurations are shown: Figure 4(a) shows the equilibrium configuration in 
which both springs are unloaded (so that they both have their natural length); 
Figure 4(b) shows an arbitrary configuration such as might occur instantaneously 
during vibrations. The displacements of the particles A and B from their 
equilibrium positions will be called respectively, x, and xg. The direction indicated 
by the single arrowheads represents the direction of displacement, velocity, 
acceleration and force which I propose to consider as positive. The choice of 
positive direction is, as usual, arbitrary, but must always be made clear and must 
be adhered to consistently in all subsequent work. 


The first thing to do is to write down the equation of motion for each particle. In 
order to do this we must know the forces due to the springs, in any arbitrary 
position such as that shown in Figure 4(b). Looking first at the spring of stiffness 
k, we find that in going from Figure 4(a) to Figure 4(b) its right-hand end has 
moved a distance x, to the right and its left-hand end has remained fixed. The 
spring has therefore been stretched by an amount x, and exerts a force of 
magnitude k,x, acting to the left on particle A. Using our sign convention we can 
sum this up by saying that the spring exerts a force —k,4x, on particle A. Now for 
the other spring. In going from Figure 4(a) to Figure 4(b) the right-hand end of 
the spring moves to the right a distance x, and its left-hand end moves to the 
right a distance x,. The spring is therefore stretched by an amount (xg — x4) and 
so exerts a force k(x 3 — x4) which acts to the left on particle B and to the right 
on particle A. Remembering that positive forces act to the right we can say that 
the total spring force on particle A is {kg(xg — x4) — k4x,4} and that on particle B 
is {—k (xg — X,)}. Since in our model the spring forces are the only forces acting 


By equilibrium position | 
mean the static position (i.e. 
the position at which there is 
no net force acting on the 
particle). 


This method of defining the 
positions of the particles is 
different from the one used in 
previous units where all 
positions were referred to a 
single origin. The convention 
used in this unit is much more 
convenient for our present 
purposes. 


0.1 kg 


0.2 kg 


Figure 3 


particle A 


particle B 


particle B 


(b) particle A 


Figure 4 


MST204 24.2 9 


on the particles in the direction of motion we can now write down the equations 
of motion of both particles by invoking Newton’s second law: 


kp(Xp — X4) — KaX4 = MgXy (1) 
—Kp(xp — X4) = MpXz. (2) 


After dividing the first equation by m,, the second by mg and then rearranging, we 
can rewrite (1) and (2) in the matrix form 


2 (k4 + kg) kg 
ma Ma|| x, X4 
} a 3 
kp _ ke ag ia ©) 
Mp Mp 


Thus (3) is a system of differential equations which describes the behaviour of the 
model in Figure 4. 


The next task is to find the solutions of (3). That is, find expressions for x, and xz 
as functions of time. There are, of course, many different solutions of (3) but in the 
next subsection we shall show that it is possible to find some particularly simple 
solutions, called normal modes, with the property that both particles execute 
simple harmonic motion at the same angular frequency. The importance of these 
normal mode solutions will become apparent at the end of this section where we 
show that any solution of (3) can be expressed as a sum of normal modes. 


2.2 Normal modes 


One way of solving the equation of motion (3) would be to turn immediately to 
the method described in the unit on systems of differential equations. Using this 
method, one can obtain the general solution of the equations of motion. However, 
before doing this I shall try to gain a more intuitive feeling for the behaviour of 
the system by considering some special solutions which can be obtained without 
reference to the simultaneous differential equations unit. To find them, let us recall 
that the typical motion of an undamped vibrating system with only one degree of 
freedom is simple harmonic motion, and ask ourselves whether a similar motion is 
possible with two degrees of freedom. More specifically, we ask ourselves the 
question: is it possible for the two-degrees-of-freedom system in Figure 1 to move 
in such a way that both particles execute simple harmonic motion with the same 
angular frequency? To answer this question let us see whether it is possible to 
satisfy Equation (3) with sinusoidal functions: 


x,4(t) = Asin(@t + 4) 
x(t) = Bsin(wt + dg) 


in which A, @4, B, dg and w are constants to be determined. Without loss of 
generality we shall impose the constraints that both A and B (the amplitudes) and 
w (the angular frequency) are positive and that both ¢@, and ¢z lie between —7z 
and 7. 


Using the above expressions for x, and xz, it follows that x, = Awcos(a@t + 4) 
and Xp = Bwcos(wt + dg) and, differentiating a second time, that 

X4 = —Aow?’sin(wt + d,) and X, = — Bw’ sin(wt + dg). We can rewrite the last 
two expressions more simply as follows: 


X4= —W°x4 (4) 
Xp == a hg: (5) 


Equations (4) and (5) are, of course, merely alternative statements of the fact that 
we are assuming that the two particles move with simple harmonic motion about 
their equilibrium positions, the frequency being the same for both particles. In 
matrix form Equations (4) and (5) can be written 3 


esl- ot} 
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Substituting this expression into Equation (3) gives 


= (kg + kp) kp 
Mma, Mma, XA So ae XA . 
.~ ea 6 
mp Mp 


From the unit on eigenvalues we know that Equation (6) has a solution provided 
—” is an eigenvalue of the matrix on the left-hand side. So, for each eigenvalue 
—q* there are solutions of Equation (3) in which both particles are oscillating 
with the same angular frequency w. Such solutions are called normal modes and 
the angular frequency of the motion is called the normal mode angular frequency. 


Finding the normal mode angular frequencies 


To avoid an excessive amount of algebra I shall illustrate how to calculate the 
normal mode frequencies by working through an example with particular values 
for m4, Mg, k, and kz. 


Example 1 

Suppose that the parameters for the model shown in Figure 1 are m, = 6, mg = 4, 
k, = 20 and kg = 10. Calculate the normal mode angular frequencies of the model. 
Solution 


Substituting the values of m4, mg, k, and kg into Equation (6) gives 


=> [Oi aXe) — | *4 
20° =—25 XB = XB 


Since the matrix is 2 x 2 its eigenvalues are easily found by solving the 
characteristic equation. Here I shall call the eigenvalues —q’, rather than / as in 
the unit on eigenvalues, and so the characteristic equation takes the form 


Se eo 1.67 
25 —25-4+ wr 


Since it is the values of the angular frequencies w that interest us, rather than the 
eigenvalues — 7, we can proceed as follows: we expand the determinant to give 
the equation 


(w?)? — 7.5(w7) + 8.33 = 0, 
and then solve for w? by the usual method for quadratics. This gives 


15+ STS 45 a 
= 2 


= 1.36 or 6.14, 


and so the possible values for the angular frequency w are 


@ = ./1.36 = 1.16 “and w= ./6.14= 248. 


There are therefore two possible frequencies at which both particles of the two- 
degrees-of-freedom model can move with simple harmonic motion. More briefly, 
the model has two normal mode angular frequencies. I shall make it a rule 
throughout this unit to call the lowest value normal mode angular frequency @, 
and to number the other values in order of magnitude. The normal mode with 
angular frequency «, is called the first normal mode, the one with angular 
frequency w, is called the second normal mode — and so on. In this case we have 


w, = 1.16 (rad s~') and W>, = 2.48 (rads~'). 


Ww? 


Now try some similar calculations for yourself: 


Exercise 3 


Assuming that in the model shown in Figure 4 the stiffnesses are k, = kg = k and the 
masses are m, = mz = m, calculate the two normal mode angular frequencies. Do the 
calculations as accurately as you can but give the numbers in the final result accurate to 
two places of decimals. 


[Solution on p. 31] 


All calculations in this unit 
are worked to the full 
accuracy of a calculator, but 
are recorded here to two. 
decimal places for 
convenience. 
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Exercise 4 


Suppose that the parameters for the model shown in Figure 1 are m, = 1.67 x 10°, 

mg = 167, k, = 12 x 10° and kg = 5.9 x 10°. Calculate the two normal mode angular 
frequencies. Do the calculations as accurately as you can but give the final results to the 
nearest integer. 


[Solution on p. 31 | 


To sum up what we have done so far: we began by investigating the conditions 
under which the two particles of the model in Figure 4 could both move with — 
simple harmonic motion at the same angular frequency. We did this by looking for 
special solutions of the equation of motion (3), having the form 


X4(t)| | Asin(wt + Pg) (7) 
x,(t)| | Bsin(at + dg) 
and found that there are two angular frequencies (w = w, and w = w>) for which 
such solutions exist. The motions described by these special solutions are called 


normal modes and the angular frequencies w, and w, are called the normal mode 
angular frequencies. 


Although we have seen how the normal mode angular frequencies w, and w». are 
found, as yet we have no information concerning the amplitudes and phases of the 
normal modes. We turn our attention to this now. 


Normal mode displacement ratios 


To a certain extent the amplitudes and phases of the normal mode oscillations are 
arbitrary. They depend on the starting positions and velocities of the particles and 
cannot be determined without additional information. What we can do however is 
find a relationship between the phases 4, and 3, and between the amplitudes A 
and B, when the particles are executing normal mode motion. To do this we go 
back to the conditions for normal mode motion given by the eigenvector equation 


(6): 


hat he) ke 
i ma | | Xa(t) | _ Soe X(t) 
ke _ke BAP ie bot (6) 
Mp Mp 


We have seen that normal mode motion occurs when — «7 is an eigenvalue of the 
matrix on the left-hand side. But so far we have not considered the eigenvectors. 
For normal mode motion the particles must move in such a way that 

[x4(t) X,(t)]’ is always equal to an eigenvector corresponding to the eigenvalue 
—w*. This means that x,(t) and x,(t) must satisfy the equations 


[ se Catto (0 = Eg) — (0) (7) 
Ma, Ma, 

£8 y(t) + [ 0 (8) 

Mp MB 


From either of these equations we obtain 
X p(t) = Rx a(t) (9) 
where R is a constant given by 


pa ket ka = mo? _ kp 
= kp kg — mg? 


The first expression for R comes from Equation (7); the second from Equation (8). 
Provided w is a normal mode angular frequency the two expressions are equal. 


From (9) we have 


x 4(t) = x 4(t) ee 1 
fon . Fadl 2 salt} R| 
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So for normal mode motion [x,(t) x ,(t)]’ is equal to the eigenvector [1 R]’ 
multiplied by the sinusoidal function x4(t) = Asin(@t + 4), or in terms of the 
component functions x,(t) and x,(t): 


x*,(t) = Asin(wt + d,4) 
— Xp(t) = RAsin(wt + dy). 


Since R is a constant, Equation (9) tells us that the two displacements x,4(t) and 
x,(t) are proportional, and that the ratio x4(t)/x,(t) (for xg(t) non-zero) is equal 
to R for all values of t. For this reason R is called the displacement ratio of the 
normal mode. 3 


Because x,(t) and x,(t) are proportional, the two particles pass through their 
static positions at the same time. There are two possible types of motion, 
depending on the sign of the displacement ratio R. If R is positive, then (9) shows 
that x, and x, always have the same sign so that the two particles are always 
both on the same side of their static positions. In this case the two particles are 
said to be in phase (see Figure 5(a)). On the other hand if R is negative, then x4 
and xp, always have opposite signs, so that both particles are always on opposite 
sides of their static positions. In this case the two particles are said to be phase- 
opposed (see Figure 5(b)). 


— — ine ih - eee Se emer CC 


Asiilial hau” 


A sin (wt + d4) 


B sin (wt + dp) =. 
' | ae eae B sin (wt + dp) 
O45 = Op | - 
(a) In-phase motion (b) Phase-opposed motion 


Figure 5 


By looking carefully at the graphs in Figure 5 we can establish a relationship 
between the amplitudes A and B and between the phases d, and @3. We start by 
considering amplitudes. If we stretch the graph of x, in each of Figures 5(a) and 
5(b) by the factor |R| (that is, by the magnitude of R) then it will have the same 
amplitude as the graph of xg. This observation enables us to write down the 
following relationship between the amplitudes: 


B = |R|A, 
or 
B 
R = pre 
RI=5 


So the magnitude of R gives the ratio of the amplitudes. 


Let us now turn our attention to the phases. It is convenient to consider the cases 
of in-phase motion and phase-opposed motion separately. If the particles are in 
phase as in Figure 5(a) then stretching the graph of x, vertically by |R| will 
actually make it coincide with the graph of xg. For this to happen both 
oscillations must have the same phase. So in this case we have 


Pa oz Pp. 


If, however, the particles are phase-opposed then, in addition to stretching the 
graph of x, vertically by |R|, we will have to shift it sideways a distance z to make 
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it coincide with the graph of xg. So we conclude that the phase of xg must differ 
from the phase of x, by z. That 1s, 


le — Pal = 7. 
In the expression x,(t) = RAsin(wt + $4) given above, the phase difference is 
automatically taken care of by the sign of R. If, however, we want an expression 
for x,(t) explicitly in terms of amplitude and phase we would have to write 

X,(t) = |R|A sin(@t + dz) 
where @g = od, if R>O or lbp — d4| = 2 (with —-nt<gdgsnj)ifR<0O. 
To summarize, each normal mode has a displacement ratio that specifies the 
relationship between the two particles: the magnitude of the displacement ratio 
gives the ratio of the amplitudes of the oscillations; the sign of the displacement 
ratio determines whether the oscillations are in phase or phase-opposed. If we 
denote the displacement ratio for the first normal mode by R, and the 
displacement ratio for the second normal mode by R, then the motion of the 
particles in the first normal mode is given by 

X4 = A, sin(w,t + 4), Xp = A,R,sin(@,t + G1) 
where A, and @, are constants (@; was previously called @,). The motion in the 
second normal mode is given by 


X4A= A, sin (wt + P2); XZ >= A,R, sin (>t =} p>). 


Finding the normal mode displacement ratio 


As an illustration we shall work out the normal mode displacement ratios for the 
model in Example 1. 


Example 2 


In Example | we considered the model shown in Figure 1 with m, = 6, mg = 4, 
k4 = 20 and kz = 10. 


Calculate the normal mode displacement ratios for this model. 
Solution 


In Example 1 we considered the eigenvalue problem 


—5 1.67 XA ee ae XA 
2.55 —2.5 XB ee XB 


and found two normal mode angular frequencies w, and w,, where w? = 1.36 and 
«3 = 6.14. To find the corresponding normal mode displacement ratios R, and R, 
we look for eigenvectors of the form [1 R,]’ and [1 R,]’. We do this in the 
usual way by substituting the appropriate values of w* into the equation 


—5+? 1.67 447,19 9) 
2.5 =—25 fay a7 78 
and solving for x, and xg. 


For the first normal mode we put w? = 1.36 into (9) to obtain the simultaneous 
equations 


—3.64x, + 1.67x, =0 
2.5x4 — 1.14x, = 0. 
Both these equations give 
Xp = 219 X, (to two places of decimals) 


and so putting x, = 1 we obtain the required eigenvector [1 2.19]’. It follows that 
the normal mode displacement ratio corresponding to the normal mode angular 
frequency w, is R, ~ 2.19. 
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For the second normal mode, we put w* = 6.14 into (9) to obtain the 
simultaneous equations 


1.14x, + 1.67x, = 0 
2.5X4 + 3.64x3 = 0. 
Both these equations give 
Xp = —0.69x A - (to two places of decimals) 


and so putting x, = 1 we obtain the required eigenvector [1 —0.69]’. It follows 
that the second normal mode displacement ratio is R, = —0.69. The negative sign 
means that the motions of the two particles are phase-opposed so that when one 
particle has reached the extreme point of its travel to the right of its static 
position, the other particle will be at the extreme point in its travel to the left of its 
static position. It also means that during this motion the velocities of the two 
particles will always be in opposite directions (except, of course, when the 
velocities are instantaneously zero at the extreme ends of the travel of each 
particle). 


The general undamped lumped-parameter model with two degrees of freedom 


In this subsection I have dealt almost exclusively with the model shown in Figure 
1 (page 7). However, all that I have said can be applied equally well to any 

undamped lumped-parameter model with two degrees of freedom. This is because 
the equations of motion for such models can always be written in the matrix form 


Hx = x, (10) 


where H is a 2 x 2 matrix, and this is a straightforward generalization of the 
equation of motion (3) for the model in Figure 1. Indeed, we can obtain (3) from 
(10) by putting 


tk te 

= M4 M4 
’ 
Mp Mp 


The techniques developed so far in this subsection have been based on Equation 
(3) but all these techniques can be applied to the more general Equation (10), as 
described in the following procedure. 


Procedure 2.2: To find the normal modes of an undamped lumped- 
parameter model with two degrees of freedom 


1. Write down the equations of motion for the system and express 
them in the matrix form 


Hx = x. 


2. Write down the characteristic equation of the matrix H in 
terms of —w? and solve’for w to obtain the two normal mode 
angular frequencies w, and @. 


3. Find corresponding eigenvectors in the form [1 R,]’ and 
[1 R,]' to obtain the two normal mode displacement ratios R, 
and R>. 


4. The motion of the particles in the first normal mode is given 
by 


X4 = A, sin(w@,t + 1) and Xp = R,A,sin(@,t + ¢;) 


and the motion of the particles in the second normal mode is 
given by 


X4 = A, sin (wt + 2) and Xp = R,A, sin (wt + P 2). 
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You may wonder what happens if the matrix H in this procedure has positive 
eigenvalues, for it would then be impossible to find real values for the normal 
mode angular frequencies. It turns out, however, that the matrices H resulting 
from lumped-parameter models can never have positive eigenvalues. The only 
complication that can arise is for one of the eigenvalues to be zero. We shall deal 
with such cases in Section 4. 


Exercise 5 


Using the data and results of Exercise 4 determine the corresponding values of the normal 
mode displacement ratios. 


[Solution on p. 32] 


Exercise 6 
Calculate the normal mode displacement ratios for the model in Exercise 3. 


[Solution on p. 32] 


2.3 The general motion 


So far, all that we have done has concerned the normal mode motion of simple 
idealized systems like the one shown in Figure 1 (on page 7). Normal mode 
motion means that both particles move with simple harmonic motion at the same 
angular frequency, though generally with different amplitudes. But suppose you 
actually made up a system like the one in Figure 1, displaced each of the two 
particles by an arbitrary amount from its equilibrium position and then released 
both particles: you would probably find that neither particle moved with simple 
harmonic motion. Indeed the motions would probably look quite irregular and 
complicated. The television programme associated with this unit includes a 
demonstration of this. So the question arises: Is there any connection between the 
kind of motion I have just described and the normal modes? The answer is ‘yes’; 
and, what is more, if we know the normal mode frequencies and displacement 
ratios, then we can predict the motions of the particles which result from any 
arbitrary set of initial conditions. 


To justify this statement we go back to the equations of motion. From the last 
subsection we know that the equations of motion for an undamped lumped- 
parameter system with two degrees of freedom can be written as a system of 
differential equations of the form 


Hx = x. (10) 


This type of system was discussed in the unit on the systems of differential 
equations, where (for negative eigenvalues) the general solution was shown to be 


X= . a.(C,cos(./ —A,t) + D,sin(./ — A, t)). 


r=1 


In this solution n is the number of equations in the system, 4, is the rth 
eigenvalue of H, a, is the rth eigenvector of H and the C,’s and D,’s are arbitrary 
constants. We can now rewrite this general solution so as to conform with the 
notation developed during our discussion of normal modes in the last subsection. 


In the last subsection we confined our attention to models with two degrees of 
freedom. For such models n = 2. Also, we used the notation —w,” for the 


eigenvalue A, and so we can replace ./—4A, by w,. Furthermore the eigenvector 
corresponding to the eigenvalue —«,? was chosen to be of the form [1 R,]* 
where R, is the rth normal mode displacement ratio. Thus in the notation of this 
unit the general solution of (10) is 


1 1 
X= & Ke cos w,t + D,sin@,t) + Fe (cs cos wt + D,sin@,t) 
1 2 


or in terms of x, and xz: 
x4(t) = (C, cos@,t + D, sinw,t) + (C,cos@2t + Dz sin@ Et) (11) 


xXp(t) = Ry(C, cosw,t + D, sinw,t) + R2(C,cos@.t + D,sinw ft). 


See Unit 22, Section 3. 
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Now we know from Section 1 that the sinusoidal functions in the brackets can all 

be rewritten to give the general solution in the alternative form: 
x,(t) = A, sin(w,t + @,) + A2sin(@2t + d2) (12) 
x p(t) = R,A, sin (@,t + d;) + R,A, sin (>t + >). 


Both (11) and (12) are useful forms for the general solution of (10). The form given 
in (12) shows clearly that the general solution is a sum of two normal modes: one 
with angular frequency w, and displacement ratio R,; the other with angular 
frequency w, and displacement ratio R,. Now the sum of two sinusoids of 
different angular frequenciés is not a sunusoid and so by choosing different 
combinations of amplitudes A,, A, and phases ¢,, @2 we can obtain an unlimited 
number of non-sinusoidal solutions for x, and xg. Thus for an arbitrary set of 
initial conditions we should, in general, expect a non-sinusoidal solution. The only 
exceptions are those particular initial conditions which satisfy the requirements of 
normal mode motion, and hence cut out the other normal mode angular frequency 
altogether. In general, however, the motion will be a non-sinusoidal combination 
of normal mode terms. 3 


If we know enough about a system to be able to find the normal mode angular 
frequencies w,, 2 and the normal mode displacement ratios R,, R2 then we have 
enough information to write down the general solution. However, if we want to 
find a particular solution we still have to find four arbitrary constants (either C,, 
D,, C2, D> in the case of (11) or A,, Az, @1, 2, in the case of (12)) and so we 

_ need four initial conditions. Here is an example. 


Example 3 


For the same data which we used in Exercises 3 and 6, i.e. kg = kg = k and 
m, = Mz = mM, predict the motion of each particle when the system is released with 
the following initial conditions 


x4(0) = X, x,(0) = 0 and x4(0) = X,(0) = 0 


where X is an arbitrary number. 


Solution 


From Exercises 3 and 6 we know that the normal mode angular frequencies are: 


k k 
oO; = 02, | and @> = 1.62 ‘ 


and the normal mode displacement ratios are: 
K, = 62 and R, = —0.62. 
We shall use the form of the general solution given in (11): 
x4(t) = (C, cos@,t + D,sinw,t) + (C,cos@ t+ D,sinw 1) ) 
Xp(t) = Ry(C, cosm,t + D, sin@,t) + R2(Czcos@t + D2 sina Lt). s ai 
We shall also need the velocities 
X4(t) = w,(—C,sinw,t + D; cosa@,t) 
+ W,(—C,sinw,t + D, cos wt) 
Xp(t) = ©,R,(—C, sina@,t + D, cosa ,t) 
+ W2R,(—C,sin@ zt + D2 cos wo). 
Substituting in the initial conditions we obtain 
C,+C,=xX 
R,C, + R,C, = 0 
w,D,+@,D,=0 
w,R,D, + @2R,D,=0 
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which can be solved for C,, C>, D,, and D,. From the second pair of equations 


we obtain (since R,; # R>) 
D, — D, SSS 0. 
From the first pair we obtain 
R,X 
R,—R2 
R,X 
mi 


0.62X 
1.62 + 0.62 


_ _1.62X 
~ 1.62 + 0.62 


C,= ~ 0.28xX 


C,= ~ 0.72X. 


Substituting these values into (11), together with the values of @,, w2, R,; and Rp, 


we obtain the required particular solution: 


x4(t) ~ 0.28X cos 


062, |) + 0.72x-cos (102, [1 


val) = 045K 00s 062, Ar} — 045 oo 62, |r), 


These solutions are shown plotted in Figure 6 assuming k/m = 1 and X = 1 (only 
part of a cycle is shown). Notice that the curves are not sinusoidal, and check that 


the required initial conditions are satisfied. 


= ee 
ae 
& ae +4 2 pt yy 4 Ht +H Ht t+H 
a a ae Ha 
s cry 5 ome t co 
HME rit A rE 
SEEN ae 
Efe ttf Peet fist 
2a . F a 


Figure 6 
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Using the same data as in Example 3 above show that when the particles are released with 
the following initial conditions, both particles will perform simple harmonic motion at the 


same frequency: 
X4(0) = X, xp(0) 


Here X is an arbitrary number. 


— R,X and 


[Solution on p. 32] 


Summary of Section 2 


X4(0) = X(0) = 


This section deals with undamped lumped-parameter models with two degrees of 
freedom. Such models have two angular frequencies at which both masses can 
move with simple harmonic motion. This type of motion is called a normal mode 
and the angular frequencies are called normal mode angular frequencies. In each 
normal mode the displacements of the two particles are proportional. The constant 
of proportionality is called the normal mode displacement ratio. A positive 
displacement ratio means that both particles move in the same direction during 
normal mode motion (in-phase); a negative displacement ratio means they move in 


opposite directions (phase-opposed). 


The general free motion of an undamped lumped-parameter model with two 
degrees of freedom (i.e. the motion which results when the particles are released 
from arbitrary initial conditions) is not sinusoidal but can be shown to consist of a 
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superposition of appropriate proportions of two normal modes. The general 
expressions for the displacements x, and x z of the two particles take the form 


X4 = A,sin(w,t + ¢;) + Az, sin(wot + g2) 

Xp = R,A,sin(@, + ¢;) + R2A2sin(wot + d2) 
where w, and w, are the normal mode angular frequencies, R, and R, are the 
corresponding normal mode displacement ratios and A,, A, @; and @2 are 


constants which satisfy A, > 0, A, > 0, —xt< ¢, <n and —7z < ¢2 < x but which 
are otherwise arbitrary. 


3 Vibration absorption (Television Section) 


3.1 Pre-broadcast notes 


The programme is about undamped lumped-parameter systems with two degrees 
of freedom and it refers to the results of a generalized analysis as follows: 


Figure 1 shows the model and the notation used in the programme. The 
displacements x, and xg are measured (as usual in this unit) from the equilibrium 
positions, so that the spring forces available to accelerate the particles (over and 
above the spring forces required to support the particles) are k,x, and kp(xp — x4) 
respectively. 


Taking the direction of these forces into account, it follows that the equations of 
motion are 


kp(xp — X4) — kaX4 = MgX4 (1) 
—Kp(xg — Xa) = MpXz. (2) 
We can find the normal modes by substituting ¥, = —w?x, into (1) and 
Xp = —W?xg into (2) to obtain 
kg(Xp — X4) — k4X4 = —M4yw°x4 
—k»(xp — Xa) = —mM,W"X>p. 


These equations can be rewritten in the form of an eigenvector equation: 


a (Kg + kp) kp 
7 Bg Ma Xa] _ 2| Xa 
= —q@ ; 
or 
mp Mp 


or alternatively, 


MU ee 
: Ma Mma, ine = Hl (3) 
Mp Ms 


The matrix on the left hand side of this equation is quoted and discussed in the 
programme. 


Exercise 1 
Show that the characteristic equation of the model in Figure 1 is 
k katk kak 
=e [ie 5 Ma = Ae) an 4 Sate 9 
Mp my, mamMp 


[Solution on p. 32] 


From the characteristic equation found in Exercise | it is possible to find the 
values of w*, each of them corresponding to one of the two normal mode angular 
frequencies w, and @. 


equilibrium 
position of 
particle A 


equilibrium 
position of 
particle B 


Figure 1 
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Now view the television programme ‘Vibration absorbers’. ° 


3.2 Vibration absorbers 
TV24 


The programme demonstrates free and forced vibration in undamped lumped- 
parameter models with two degrees of freedom, and-discusses the way in which 
these models can predict the performance of real systems, and help with the design 
of vibration absorbers whose purpose is to suppress unwanted vibrations. 


To illustrate free vibrations, there is a studio demonstration of the fact that, given 
the right initial conditions, it is possible for the free motion of a system like that in 
Figure 1 to take the form of simple harmonic motion of both particles, the 
frequency being the same for both. But it is also shown that, in general, for 
arbitrary initial conditions, the free motions of the particles are not simple 
harmonic. 


Exercise 2 


The studio demonstrations use a double spring-mass system similar to that shown in 
Figure 1 in which 


Ma = 0.38(kg) ka = 15(Nm~*) 
mz = 0.146(kg) kp = 5.8(Nm~'*). 


Estimate the normal mode angular frequencies for this system assuming it to be undamped y Uy 
and assuming the springs to be perfect. os ass 


[Solution on p. 33] 


Exercise 3 

For the demonstration system specified in Exercise 2, estimate a ratio for the initial 
displacements of the particles from their equilibrium positions so that, after release, both 
particles move with simple harmonic motion at the same angular frequency. 


equilibrium positions 
of the particles 
when there is 

no forced motion 


[Solution on p. 33] 


In the programme, forced motion is illustrated by giving the upper end of the 
spring in Figure 1 a vertical sinusoidal motion y = YsinQt. In other words, what 
was the fixed end of the system is now made to move with simple harmonic 
motion. The studio demonstration shows that after a short time both particles 
settle down to a steady-state consisting of simple harmonic motion with the same 
angular frequency as the motion imposed on the upper end, i.e. with angular 
frequency 9. The equations of motion can be derived from Figure 2 (the Figure 2 
displacements x, and xg of the particles being measured from the same positions 

as in Figure 1): 


—K4(x4 — y) + kp(xp — X4) = myXy : (4) 
—kp(xg — X4) = mpXz. (5) 


Notice that, as pointed out in the programme, these two equations are very similar 
to Equations (1) and (2) which were obtained for the free motion of the system. 
The only difference is the presence of the y in Equation (4). 


Remembering that y = Ysin 1, we can rewrite the last two equations as follows: 
MaX4 + (kg + kp)x4 — kpxp =k, YsinQt (6) 
MpXp+kpxp—kpxy a (7) 

We can find the steady-state solution of these simultaneous equations by following 


the procedure explained in Subsection 4.2 of Unit 22. 


Exercise 4 
Find the steady-state solution of Equations (6) and (7). 


[Solution on p. 33] 
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From the solution to the last exercise it is clear that the amplitudes (or maximum 
values) of x4 and xg, are, respectively, 


A = (kp a“ mpQ’)ksY 
(ky + ky — mgQ?)(kg — mpQ?) — kp” 
and 
= kakpY 
(ky + kp — m4Q)(kg — mpQ?) — kg” 


Figure 3 shows these amplitudes plotted against 2. 


B 


amplitude of x4 amplitude of xz 


0) w > 9) 0 w w> 
Figure 3 (a) main mass (b) absorber mass 


The two noticeable features of the graphs are: 
(i) that there is a value of Q for which x4 = 0; 


(ii) that there are two values of Q near which both amplitudes become very large. 
These values of Q are known as the resonant angular frequencies. 


Resonance will occur at the values of Q for which the denominator of the 
expressions for the amplitudes A and B is zero. (Note that both expressions have 
the same denominator.) The condition for resonance therefore is 


(ka + Ky = m4Q? (kg — mgQ?) — kz? = 0, 


1.€. k ake ae k pm ,Q? me k 4mpQ? me k pm pQ? + m,4mpQ* = 0, 
ce (" ‘thal gp 4 Tae 0 
Mp ma, m,mMpg 


If we put Q equal to the normal mode angular frequency «, we get the 
characteristic equation of the system which you derived in Exercise 1. It follows 
that the normal mode angular frequencies are also the resonant angular 
frequencies of the system. 


The value of Q for which x, = 0 can clearly be found by putting the numerator of 
the expression for A equal to 0. Since k,Y 4 0 (otherwise we should have no 
forced motion), we must have 


kp are mpQ? = 0, 
k 
i.e OQ? = —*, 
Mp 
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At this angular frequency particle A is stationary even though the top end of the 
system continues to move according to y = Ysin Qt, and particle B continues to 
vibrate, too. Under these circumstances the system is said to act as a vibration 
absorber, m, being the main mass whose vibrations are absorbed (so that it is 
stationary) and mg, being referred to as the absorber mass. As you can see from the 
graph, fine tuning is essential for good vibration absorption; even quite a small 
change in the value of Q can move the value of x, quite a long way away from 
zero. 


The vibration absorber designed by Southampton University to reduce the 
vibrations of the car deck on a Sealink ferry (now called Earl William, formerly 
Viking II) works on the principles I have just outlined. The vibrations are due to 
the action of the propellers. The mass and the stiffness of the car deck constitute 
the main mass my, and the stiffness k,. A mass of water supported on a column of 
air makes up the absorber mass and its spring. A schematic sketch of the 
arrangement is shown in Figure 4. 


In the programme, Dr Hunt of Southampton University demonstrates a small 
version of this device and points out that it might be possible in future for the 
tuning to be done automatically so as to maintain a low amplitude of vibration of 
the car deck at all engine speeds. The present design does not have this facility but 
is tuned to the vibration frequency at normal speed, which is about 13 hertz (i.e. 
13 cycles per second). 


Summary of Section 3 


The free motion of an undamped lumped-parameter model with two degrees of 
freedom can be expressed as a superposition of the two normal modes of 
vibration. 


For the forced motion of the undamped lumped-parameter model with two 
degrees of freedom shown in Figure 2 the normal mode angular frequencies are 
also the resonant angular frequencies. The graph of steady-state amplitude against 
angular frequency for each particle therefore shows two resonant angular 
frequencies. Between these there is one angular frequency at which the amplitude 
of one particle (the one with mass m, in Figure 2) is zero though the vibration of 
the other particle has a finite amplitude. This angular frequency is given by 


[k 

Q= |/—* where kg is the stiffness of the spring connecting the two particles. At 
Mp 

this angular frequency the system is said to act as a vibration absorber. This 

principle is used to suppress unwanted vibrations, as in the case of the car deck on 


the ship featured in the television programme. 


4 Extending the scope 


4.1 About this section 


In this section, I return to the discussion of normal mode motion. So far in this 
unit I have concentrated mainly on the simplest model which demonstrates the 
properties of normal modes: the undamped, two-degrees-of-freedom, lumped- 
parameter model shown in Figure 1 of Section 2. In the first part of this section I 
want to consider a physical system which is similar to this model (not least 
because it, too, has two degrees of freedom) but which takes a slightly different 
form. 


The rest of this section will be devoted to coping with more than two degrees of 
freedom. For this, some of the mathematical methods from the unit on eigenvalues 
will again be very useful. I shall again stick to undamped models and I may as 
well say here and now that those with more than two degrees of freedom behave 
in much the same general way as those with just two, except that they have more 
normal mode angular frequencies and normal mode displacement ratios and 
require bigger matrices which are somewhat more laborious to handle. But that is 


Figure 4 


mass of water 


stiffness of 
air column 


mass of car deck 


stiffness of 


car deck 


forcing due 
to propellers 
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merely a matter of technique—the treatment of more than two degrees of freedom 
is really just an extension of methods we have already used: no new principles are 
involved. 


4.2 A degenerate model 


The model I want to discuss next is shown in Figure 1. You have already met this - : 
model in the unit on simple harmonic motion — the first unit on vibrations — but 

you may find it interesting to meet it again in a slightly different context. It is 
different from the two-degrees-of-freedom models that I have dealt with so far in 
this unit, in that no part of the model is attached to a fixed point. 


Figure | 


We shall treat the model in the same way that we have treated the other models in 
this unit. We begin by choosing an equilibrium configuration from which to 
measure the displacements of the particles. Because no point of the model is fixed, 
there are many ways of doing this. Any configuration in which the spring has its 
natural length /o will do, for there will then be no forces acting on the particles. It 
does not matter which equilibrium configuration we choose so long as we pick a 
particular one and stick to it. Figure 2 shows a snapshot of the model at some XA XR 
instant after the particles have been set in motion, together with the equilibrium ee \ | 
positions from which the displacements x, and xg are to be measured. The | | 
arrowheads for x4 and xg define, as before, the positive direction for 
displacements, velocities, accelerations and forces. The spring force is k(xg — x,), 

and so, taking the direction of forces into account, the equations of motion of the Figure 2 
two particles are: 


k(xp — X4) = M4X4 
—k(xp — X4) = mgxp 


which can be written in the matrix form 


k k 
mM, mM, XA = X4 
k kode 
Mp Mp 
For normal mode motion we put ¥, = —w?x, and Xp = —w7’x, to obtain the 
eigenvector equation 
k , k 
M,° Ma, in -|5| (1) 
k [ k aaa Ww XB 0 
Mp Mp 


The characteristic equation of the model is therefore 


k k k? 
m 


mM, B m,mMpg 
2 2 
4 2 2 k k = 
or w™ — —w* — —@* + —— —- = 0 
mM, Mg msm, M,sMpZ 


: k 
i.e. ao ie + Jo" = 0. 


This is fairly easy to solve because there is no constant term and so we can factor 
out w? to obtain 
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A system for which one of the normal mode frequencies is zero is called a 
degenerate system. The effect that this has on the normal mode motion can easily 
be seen by substituting w = 0 into one of the equations in (1) above. From the top 
equation we get 


k 
— —x, + —x, = 0 
Ma Ma 


SO Ky Xp. 
Thus, at any instant of the normal mode motion, both ends of the spring are 
displaced by the same amount and so this corresponds to a mode in which the 
spring is neither compressed nor extended. This means that the force on each 
particle is always zero and so the particles move at a constant velocity as though 
they were part of the same rigid body. For this reason, this normal mode 1s 
referred to as a rigid body mode and its existence is typical of a degenerate system. 
The other mode is unremarkable; you might like to work out its displacement 
ratio for yourself: 


Exercise | 
Derive the second normal mode displacement ratio of the degenerate system in Figure 2. 


[Solution on p. 33 | 


The general expression for the displacement of each particle will be similar to that 
which we obtained in Section 2, but with one important difference: since, as we 
have seen, the motion corresponding to the lower normal mode is one of constant 
velocity for both particles (the same velocity for both) it follows that for this 
normal mode 


X4 = Xp = C, 
and hence 
x, =X, = Ct + D: 
where C and D are arbitrary constants. 


The motion corresponding to the second normal mode is unchanged, so that, for 
the general motion 

X4>= Ct+D+ A, sin(@at 38 p>) 

Xp = Ct+D+R,A,Sin(@2t + 2). 


4.3 Three degrees of freedom 


As I have already said, a model with three degrees of freedom behaves in the same 
general way as a model with two degrees of freedom except that it has three 
normal mode angular frequencies and three normal mode displacement ratios. In 
this subsection I will consider two models with three degrees of freedom: first a 
degenerate model and then a non-degenerate model. 


A degenerate model 


This time, I will start with a numerical example. Consider the three particles and 
two springs shown in Figure 3. 


20 10 


XA XB Xo 
x N Y 


Figure 3 
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The masses are given in kg and the spring stiffnesses in Nm‘, as usual. The 
displacements are measured from the position which each mass takes up when the 
whole system is in one of its equilibrium configurations (that is, a configuration in 
which both springs have their natural lengths). We now imagine the three particles 
to be set in motion, and we proceed in the usual way. First, we need the equations 
of motion for the three particles. 


Exercise 2 
Write down the equations of motion for the three particles in Figure 3. 


[Solution on p. 33 | 


You can now write down the eigenvector equation for normal mode motion, by 
incorporating the conditions X4 = —w?x,4, Xp = —@7xg and X¢ = —w*Xxc¢ which 
specify the simple harmonic motion of the three particles. 

From this eigenvector equation you can derive the characteristic equation for the 
model in the same way as we did for models with two degrees of freedom. 


Exercise 3 


Write down the eigenvector equation for normal motion and hence derive the characteristic 
equation of the model in Figure 3. 


[Solution on p. 33 | 


So the characteristic equation of the model is 
o° — 11.5o* + 25er7 =0 


where w represents the normal mode angular frequency in rads ~ ' This is a cubic 
equation in w’, but notice again that the constant terms have cancelled, so that we 
can write 


w?(w* — 11.5m? + 25) = 0. 
Hence, either wm” = 0 or 
(wo* — 11.50? + 25) = 0. 


So we need only to solve a quadratic in w? to find the other two normal mode 
angular frequencies. We have 


7 te 
Hence 
wo," =0 if - We) 
>” = 291 i.e. w> = 1.71 (rads~') 
3" = 8.59 i.e. Ww; ~ 2.93 (rads~'). 


The lowest normal mode angular frequency is zero, and so the model is 
degenerate. We can quickly check that the lowest normal mode is a rigid body 
mode as follows. 


From Exercise 3 we know that the equations for the normal mode eigenvalue 
problem can be written 


(-4+07)x,+ 4x, =0 
3.33x, + (—5 + w*)xg + 1.67x- =0 
99 (25 ee 0. 


For the normal mode frequency w, = 0 it is clear from the first equation that 
X4 = Xp, and from the last equation that xz = x-. So for this normal mode the 
displacements of the particles are equal at all times during the motion. This ts 
therefore a rigid body mode. Since both springs remain at their natural lengths, 
the force on each particle is zero and so all the particles move with the same 
constant velocity. 
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It is interesting to characterize this rigid body motion in terms of displacement 
ratios. When we were dealing with models having two particles, each normal mode 
had a displacement ratio x,/x,4 associated with it which related the motion of 
particle B to the motion of particle A. Now that we are dealing with a model 
having three particles, each normal mode requires an extra displacement ratio 
X¢/X4 to relate the motion of the third particle C to the motion of the other . 
particles. In the case of the rigid body mode just considered x, = Xp = X; and so 
both the displacement ratios xp/x, and x-/x,4 must be 1. 


To complete the calculations, let us find the displacement ratios for the other 
normal modes. For the second normal mode we have w* ~ 2.91 and so the 
equations become 


—1.09x,+ 4x, ==) 
3.33x4 — 2.09x, + 1.67x- = 0 
2.5xp + 0.41x;- = 0. 


These can be solved by Gaussian elimination as in the unit on eigenvalues; but as 
in our previous eigenvector calculations in this unit it is quicker to start with the 
first equation, which gives x, in terms of x4: 
1.09 
XB = es == O72 1X5 
4 
Then we can use the third equation to express x; in terms of xg and thence in 
terms of x,: 


Re: Zo 
DWN Se Gh Prager © 
041 Xp 041 ~ 027%, 66x 4 


XC = 


The middle equation cari then be used as a check: 
3.33x4 — 2.09 x 0.27x,4 + 1.67 x (—1.66)x, ~ 0. 
So, if we take x, = 1, we obtain the eigenvector 


1 
Cnet 
= 166 


Because the entry for x, is unity, the entries 0.27 and — 1.66 give the two 
displacement ratios x,/x,4 and x-/x, for the second normal mode. The change of 
sign shows that at any instant during the motion particles A and B are moving in 
the same direction (in-phase) while particle C is moving in the opposite direction 
(phase-opposed) to A and B. 


Repeating the procedure with w? = w3”* = 8.59 we obtain 
Xp = —1.15X, 
Xc = —0.41x_, = 0.47x4. 

If we take x, = 1 we obtain the eigenvector 


1 
= 1.15 
0.47 


So, in the third normal mode the displacement ratios are —1.15 and 0.47. Hence 
particle B is the ‘odd man out’, particles A and C moving in phase with each 
other. 


26 MST204 24.4 


A non-degenerate model 


The model I shall use here is derived from the degenerate model we have just 

' considered by adding an extra spring and using it to anchor the model to a fixed 
foundation, thus removing the degeneracy. The resulting model is shown in 
Figure 4. — 


et N XN 


Figure 4 


The addition of the third spring is reflected in a change in the equation of motion 
for particle A which now becomes: 3 


— 10x, + 20(xp = x) = 5X4; 


the other equations of motion remaining unchanged. The eigenvector equation for 
normal mode motion of the model in Figure 4 is therefore 


—6 4 0 Ke = 
3.33 —S 1.67 XB = —w? Xpi, (2) 
0 2.55 —2.5 X¢ X¢ 


which is the same as the eigenvector equation obtained in the solution to Exercise 
3, with the exception of the term —6 which takes the place of — 4. By comparison 
with the solution to Exercise 3 the characteristic equation for the model is 


(—6 + w7V8.33 — 7.50? +o) — 1953 $25 +o) = 0 
or w° — 13.5m* + 53.33w? — 50 + 33.33 — 13.33m? = 0, 
i.e. w® — 13.5w* + 40m? — 16.67 = 0. (3) 


This equation has a constant term in it and we cannot therefore reduce it to a 
lower order by factoring out w* as we did before. So we are faced with the 
problem of solving Equation (3)—a cubic in w?. One way to proceed is to use 
one of the numerical methods from the earlier units of the course. To do this it is 
convenient to work with the eigenvalues 4 = —q@* than directly with the angular 
frequencies w. In terms of 4 the characteristic equation (3) becomes 


| —A3 — 13.51? — 404 — 16.67 = 0, 
or A> + 13.517 + 401 + 16.67 = 0. 


This equation was solved numerically in Unit 18 (Section 2, Exercise 1) using the 
Newton-Raphson method, where the roots were found to be 


A, = —0.49707, A, =~ —3.5464, A3 =~ —9.4565. 


We could, of course, have calculated these eigenvalues using one of the numerical 
procedures from the unit on eigenvalues, and that would have been an equally 
reasonable way to proceed. 


Once the eigenvalues have been found we can calculate the normal mode angular 
frequencies using 4 = —w*. We obtain (to two places of decimals) 


@,7=050 ie. w, = 0.71 (rads~') 
WwW,” = 3.55 i.e. W> = 1.88 (rads~') 
w3” = 9.46 i.e. w3 ~ 3.08 (rads '). 


The corresponding eigenvectors (and hence the displacement ratios) can be found 
by substituting these values of w in the eigenvector equation, as before. 
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The general case of free motion, for three degrees of freedom 


Lastly, we must look at the general motion of the particles in terms of the normal 
modes. We first consider the non-degenerate case. Just as for the two-degrees-of- 
freedom model we wrote: 


e — Aa] p, [sino + @,)+ Aal [sins + do), 


where 
@, and @, are the normal mode angular frequencies, 
R, and R, are the normal mode displacement ratios, 


A,, Az, @; and @, are constants whose values can be determined by the 
initial conditions, 


so for a non-degenerate, three-degrees-of-freedom model we can write the general 
solution in the form 


Ny 1 1 
XB = Ay Rig sin(w,t + @,) + Az R>p sin(@ zt + 2) 
Xo Ric Ro¢ 
1 
A; R3p sin (a3t + 3) (4) 
R3¢ 
where 
@,, @2 and @3 are the normal mode angular frequencies, 
oe 
R,p is the displacement ratio —* corresponding to @, (r = 1, 2,3), 
XA 
: oe ; 
R,c is the displacement ratio = corresponding to w, (r = 1,2, 3). 
A 
Exercise 4 


Use Equation (4) to write down an expression for the general free motion of the three 
particles of the model in Figure 4. 


[Solution on p. 34] 
When dealing with degenerate models the expression for the general motion given 
in Equation (4) must be modified. The way this is done is similar to the treatment 


of the two-degrees-of-freedom degenerate model in Subsection 4.2. The following 
example demonstrates this. 


Example | 
Write down an expression for the general motion of the three particles of the 


degenerate model analysed at the beginning of this subsection. 
Solution 


For this model wm, = 0, w2 = 1.71 and w, = 2.93 and the corresponding 
eigenvectors are 


1 1 L 1 1 1 1 

eo ee ee ee see Ro = ES 
Ric il r= p= a6 R3¢ 0.47 

As in the case of the degenerate two-degrees-of-freedom model in Subsection 4.2, 

the motion corresponding to the lowest normal mode (for which w, = 0) is not 


sinusoidal. As I have already explained, the motion is one of constant velocity (or 
zero acceleration). Hence, for each of the particles, 


x = Q, 


1.e. ia C, 


As for the two-degrees-of- 
freedom case, this equation 
can be justified from the work 
in the unit on systems of 
differential equations. 
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SO x= Ct+ D, 


where C and D are constants and t stands for time. The other two normal modes 
are treated just as before and so the general motion is 


XA 1 1 
Xp} = |1)}(C+Dt)+ | 0.27 | A,sin(1.71t + 2) 
Xc 1 — 1.66 
1 
+ |—1.15 | A3sin(2.93t + ¢3). 
0.47 


4.4 Four and more degrees of freedom 


By analogy with two- and three-degrees-of-freedom systems, we can see that the 
general motion of a four-degrees-of-freedom spring-mass system of the kind we 
have been considering will have the form 


X4 4 
4 
x R 
Bl = 5. 'B | A.sin(@,t + ¢,), Degenerate models are dealt 
R ° ° 
X¢ pei rc with in the usual way by 
Xp a pee A, sin(@,t + $;) by 
aa. 


where A, and @, are arbitrary constants and [1 R,g R,c R,p]' is an eigenvector 
corresponding to the eigenvalue —w,*. Because the first component of this 
eigenvector is unity the other components give the displacement ratios of the 
particles in the rth normal mode. The eigenvectors [1 R,g R,c R,p |’ and 
eigenvalues 1 = —w,7 (r = 1,2,3,4) can be found by one of the numerical methods 
recommended in Unit 21. A four-degrees-of-freedom spring-mass system forms the 
subject of Problem 2 in Section 5. 


The generalization to more than four degrees of freedom is straightforward. 


Summary of Section 4 


The free motion of an undamped lumped-parameter model with more than two 
degrees of freedom displays the same general characteristics as shown by the two- 
degrees-of-freedom model in Section 2 except that there are more normal modes. 
A three-degrees-of-freedom model has three normal modes, a four-degrees-of- 
freedom model has four normal modes, and so on. The general free displacements 
X4, Xp and xc of the particles of a non-degenerate three-degrees-of-freedom model 
for example, may be written: 


xy 1 1 
Xp | =A, | Rip] sin(@yt + 1) + Az | Rog} sin(@2t + o2) 
XC Ric Ry¢ 
] 
+ A31|R3, | sin(w3t + ¢3), 
R3¢ 


where 


W,,@2,W3 are the normal mode angular frequencies, 


: eS ; 
R,g 1s the displacement ratio = corresponding to w,, 
x4 


: ; See : 
R,c 1s the displacement ratio a corresponding to ,, 
XA 
A,, Az, A3, >, h2, o3 are arbitrary constants. 


Degenerate models (whose lowest normal mode angular frequency is zero) may 
have any number of degrees of freedom. 
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For a degenerate model with three degrees of freedom R,g = Ric = 1 and the 
motion corresponding to the first normal mode is one of constant velocity for all 
the three particles. The first term on the right-hand side of the equation for the 
general motion of the three particles is therefore 


b 4c + D2). 
1 


The other terms remain unchanged. 


5 End of unit problems 


Problem 1 


Figure 1 shows two particles fixed to a length of string. The string is stretched, so that it is 
in tension, and, in this condition, is fixed between two rigid supports, one above the other, 
so that the string is vertical. Figure 1 shows the equilibrium position of the system. 


By making a number of simplifying assumptions (one of which is that gravity may be 
neglected compared with the tension in the string), small horizontal vibrations of the two 
particles may be modelled as a two-degrees-of-freedom system using the following equations 
of motion: 


z fs : 
Ny eg Xgl ee 


T T ; 
7a = Xp oh: = Mpxg 


where T is the tension in the string, and x, and xg are the horizontal displacements of the 
particles from their equilibrium positions. 


For the case when m, = mz = 0.05kg and | = 0.2m, determine 


(a) the value of T required to make the lower normal mode angular frequency equal to 
300 rad s~?, 


(b) the corresponding value of the higher normal mode angular frequency, 


(c) the corresponding normal mode displacement ratios. 


For each normal mode sketch the configuration of the system when the masses are furthest 
from their equilibrium positions. 


Predict the displacement of each mass at t = 2, where t is time in seconds, if the following 
conditions hold at t = 0: 


xz = 0, Ng et and i, = Xe = 0: 
[Solution on p. 34 | 
Problem 2 
Estimate all the normal mode angular frequencies of the model shown in Figure 2, in terms 


oe 
of the ratio —, given that: 
m 


(a) the characteristic equation is 
A(A3 + 6A? + 101 + 4) =0 
where 4 = —w?m/k; 


(b) one of the roots of the characteristic equation is —2. 


Figure 2 


Derive the displacement ratios corresponding to the highest normal mode angular 
frequency. 


[Solution on p. 35 | 


Please refer to the study guide 
at the beginning of the unit. 


Figure 1 
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Problem 3 


An example of the use of an undamped absorber is furnished by the motor/generator set 
shown in Figure 3. It is used to compensate for any slight unbalance of the shaft which 
would set up a periodic force and cause vibration. Notice that the total absorber mass is 
made up of two equal parts, each attached by a leaf spring at opposite ends of the machine. 
The two leaf springs have equal stiffnesses. 


: ed 
ae absorber , ez : 


a = SeNtings 
=o bsorb 
~~ absorber colon ees 
ceca masses 


Figure 3 


This motor/generator set with its absorber was the source of the data on which Exercise 4 
in Section 2 was based. That exercise was concerned only with the vertical vibrations of the 
machine and so is this problem. 


In that exercise the main mass was 1.67 x 10° kg and it was attached to the foundation by 
springs of total stiffness 12 x 10° Nm~'. Each absorber mass was 83.5 kg and the stiffness 
of the spring by which it was attached to the main mass was 2.95 x 10° Nm '. The main 
mass is subject to a sinusoidally varying force which is due to the unbalanced shaft. 


(i) At what angular frequency will the main mass remain at rest? 


(ii) During steady running, with the absorber keeping the main mass at rest, the total up 
and down movement of each absorber mass is observed to be 0.002 m. Calculate the 
maximum vertical force on the main mass due to the unbalance of the shaft. 


[Solution on p. 36 | 
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Appendix 1: Solutions to the exercises 


Solutions to the exercises in Section 1 


I. Since r= 0, we can apply Newton’s second law to Figure 
1 as if the dashpot were not present. If x = x9 when the 
spring is unstretched, then the equation of motion is 


mx = —k(x — Xo) 
or mx+kx=kxo. 


From Unit 7 we can recognize this as the equation for simple 
harmonic motion, the angular frequency being given by 


fi 
OM = — 
m 


The angular frequency (often called the ‘natural’ angular 
frequency) is therefore 4.47 x 10° radians per second. 


2 x 10° 


~ 4.47 x 10°. 


2. The figures below show two more ways of specifying the 
configuration of the system. 


Since all these representations use two coordinates the 
system clearly has two degrees of freedom. 


3. The model in Figure 3 has seven degrees of freedom. 


4. (a), (c) and (d) all have two degrees of freedom. (b) has 
one degree of freedom. (e) and (f) both have three degrees of 
freedom. (g) has four degrees of freedom. 


Notice that the number of degrees of freedom does not 
depend on the number of forces acting on each particle—it is 
a purely geometrical notion: the smallest number of 
displacements from any fixed datum (arbitrarily specified) 
that is required to describe the configuration of the assembly 
at any moment during the motion (irrespective of what forces 
may be acting). For example, the two rigidly connected 
particles in (b) need only one coordinate to specify their joint 
position (in fact, they are equivalent to just one particle); 
and so this model has only one degree of freedom. 


5. The particle has 3 degrees of freedom. If you look back to 
what I said about Figure 4 in Subsection 1.3 you will realize 
that we now need an extra coordinate to allow for motion 
into and out of the paper. 


Solutions to the exercises in Section 2 


1. The equilibrium position corresponds to zero spring force, 
so that both springs have their natural lengths. It follows 
that the 0.1 kg mass is 0.6m from the datum line and the 

0.2 kg mass is 1.2m from the datum line. 


2. Here equilibrium must cater for the support of both 
particles so that the force in the lower spring (attached to the 
0.2kg mass) must equal 0.2 x 9.81 ~ 1.96N and the force in 
the upper spring must equal (0.2 + 0.1) x 9.81 ~ 2.94N. 
The extension (beyond the natural length) of the upper 

pos 2.94 
spring is, therefore, == 0.039 m so that the upper particle 


is 0.039 + 0.6 = 0.639 m below the support. Similarly the 


; = ee 
extension of the lower spring is a5 ~ 0.026m so that the 


lower particle is 0.639 + 0.6 + 0.026 = 1.265m below the 
support. 


3. The equations of motion are: 
k(xp — X4) — kx4 = mX4 
—K(xp — X4) = mXz 


and so the eigenvector equation for normal mode motion is 


2k k 

m m XA ee. XA ~ 
k k |= . Fa 

m m 


The characteristic equation for the model is therefore 


a k 
a m a 
o Be ee 
m 
2k k k? 
i.e. fo - 3 fur -*] © =0 
m m| om 
Be 
or He Sea 
m m 


Solving this equation for w* we obtain 
3k A 9k? 4k? 
a S 5\e 
ee 2 ee. 
2 2 am 


Denoting the lower of the two values of w by w, and the 
higher by w,, we have 


k k 
w,? = [ =e = 0,38 


see ae m 


k 5 
if. a, ~062, | (ads = 


i.e. on 162 [* ads”) 
m 


4. The equations of motion are: 
59 3 10°Og 54) 12 5 10°X) = 167 « 10° X,, 
—5.9 x 10°(x, =< X4) = 167X.p, 


and so the eigenvector equation for normal mode motion is 


179x 10° 59x 10° 
pia? cur 18) 10 1 i is I x; 
59x10 59 x 10° be Se ta 
167 ee | 


The characteristic equation for the model is therefore 


fay = 10:72 « 10°) 353 x 10° 
Bee he 10° (o? —.35,33 « 10°) 


be et = 10.72 10°)? 95.93 x 10) — 104.82 « 10° = 0. 
or w* — 46.05w? x 10° + 253.86 x 10° = 0. 


== 0): 
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Solving this equation for w* we obtain 


ip? — 46.05 x 10? + v/2120.41 x 10% — 1015.45 x 10° 
| 2 
= 23.02 x 10° + 16.62 x 10°. 
So we have 
w,* = 6.40 x 10° i.e. 
w,”? = 39.64x 10° ie. 


w, = 80 (rads *) 
W, = 199 (rads *). 


5. In Exercise 4 we considered the eigenvector equation 
— 10.72 x 10° cies ae 62 ae _ | * 
35.33. x 10° 35.32 @ 1 sk, | Xp 

and found two normal mode angular frequencies w, and @>, 


where w,? = 6.40 x 10° and w,” = 39.64 x 10°. The 
corresponding eigenvectors can be found by substituting 


these values of w into 
= 10.72 x 10° + w’ 3.53 x 10° fal = H 
- 35.33 x 10° — 35.33 x 10° +.’ | | Xp 0 
and solving for x4 and xz. 
In the case w? = 6.40 x 10° we obtain the equations 
—432 xt x3555 x 10x, =0 
35.33 x 10°x, — 28.93 « 10°x, =0 


both of which give x, ~ 1.22x,. Putting x, = 1 gives the 
eigenvector [1 1.22]’, so the first normal mode displacement 
ratio is 1.22. 


In the case w? = 39.64 x 10° we obtain the equations 
28.92 x 107*x, 4353 © 10x, =0 
35.33 x 10°x, + 4.31 x 10°x, =0 


both of which give xg ~ —8.2x,. Putting x, = 1 gives the 
eigenvector [1 —8.2]’, so the second normal mode 
displacement ratio is — 8.2. 


6. In Exercise 3 we considered the eigenvector equation 


fee 
m m XA eee XA 
k k * a eS 
m m 


and found two normal mode frequencies w, and w 2, where 
k k 

w@,? = 0.38— and w,” = 2.62—. The corresponding 
m m 


eigenvectors can be found by substituting these values of w? 
into 


and solving for x4 and xg. 


k 
In the case w* = wate we obtain the equations 
k k 
ae 1.62—x4 + —Xz = 0 
m m 


k k 
—x, — 0.62—x,z = 0 
m m 
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both of which give x, = 1.62x,. Putting x4 = 1 gives the 
eigenvector [1 1.62]", so the first normal mode displacement 
ratio is 1.62. 


k 
In the case w* = 2.62 — we obtain the equations 
m. 
k k 
0.62—x4 + —Xz = 0 
m m 


k k 
Set 1.62—Xp, a= G 
m m 


both of which give xg = —0.62x,. Putting x, = 1 gives the 
eigenvector [1 —0.62]’ so the second normal mode 
displacement ratio is —0.62. 


7. We can use the equations from Example 3: 
X4 = C,cos@,t + D,sina,t 
+ C,cos@zt + D2 sinwyt $1) 
Xp = R,(C, cos@,t + D,; sina,t) 
+ R,(C,cos@,t + D,sin@21) 
X4 = @,(—C, sin@,t + D, COS Wt) 
+ w,(—C,sinw,t + D, cosa t) 
Xp = R,o,(—C, sinw,t + D; cos@,t) 
+ R,w,(—C,sin wt + D, cos @ ft). 
Substituting the initial conditions we have: 
C,+C,=X 
R,C, + R,C, = R,X 
D,@, + D.@, = 90 
R,D,q@, + R2,D2,@, = 0. 


From the second pair of equations we obtain, as in 
Example 3, 


Dy => De =. 
From the first pair we get 
ae Ss Co = 0. 
Hence, substituting in (S1) we find: 


k 
a2 |X 
m 
k 
Xp= R,C, COS @,t = R,X cos lose | ; 


This is the first normal mode motion for the system and 
both particles perform simple harmonic motion at the same 
frequency. 


xX, = CC, cost = X cos 


Solutions to the exercises in Section 3 


1. The characteristic equation is obtained by equating the 
determinant of the matrix on the left hand side of Equation 
(3) to zero, as follows: 


kat+ke ‘ kp 
ae Oe 

A A 

=0 
af gs 
Mp Mg 
That 1s, 
kathy slike 3 kg” 6 
mM, Mp m,mMp, 2 
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kgkp 


k kgtk 
or w*— [ts 4 at Mela? 
MMs 


Mp ma, 
as required. 
2. Using the characteristic equation found in Exercise 1 and 


substituting the given values for masses and spring stiffnesses, 
we get 


Si eole Mregeuten, Or eee 
0.146 0.38 0.38 x 0.146 


ie. w* — 94.46w? + 1568.13 = 0. 


Ps 


Solving this for w? gives 


94.46 + ./8923.23 — 6272.53 
wo? = 2 _ = 47.23 + 25.74. 


2 
Hence 
w,* = 21.49 Le. 
2" = 72.97 ie. 


Ww, = 4.64 (rads~*) 
W, = 8.54(rads~'). 


3. If the initial displacements from the equilibrium 
configuration are in the same ratio as one of the normal 
mode displacement ratios the system will have normal mode 
motion, as required. So there are two possible answers — 
one for each normal mode. 


The displacement ratios can be calculated from either of the 
equations which make up the matrix equation (3). The 
bottom equation is 


which on substituting the values of k4, kg, m4 and mz from 
Exercise 2 becomes 


5.8 5.8 
mee Pree, + rag -° }xs=0, 
ie. —39.73x,4 + (39.73 — w*)xp =0. 
Hence 
xy 83 


x, Bo 


This can take two possible values, one corresponding to each 
value of w?. From Exercise 2 we have w,7 = 21.49 and 

Ww? = 72.97 so for the lower normal mode angular frequency 
the required ratio for the displacements is 


ty BF 
%;, MES oe se 


and for the higher normal mode frequency it is 


<p 39.73 cS ae as 
%, 39973 - PRS 


A positive value indicates in-phase free motion of the 
particles, a negative value phase-opposed free motion. This 
tallies with the observation of the behaviour of the studio 
demonstration system. 


4. Following the procedure described in Subsection 4.2 of 
Unit 22 we look for a solution of the form 


Ss ]=ne( Se} 


which leads to the simultaneous linear equations 


—m,Q7a, + (k,4 + kp)a, = k pap = —k,Yi (S2) 
—m,Q7az “}- k pap —kpa, a4 (S3) 
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From (S3) we have 
kp : 
azg= | —* aan (S4) 
Substituting in (S2) we have 


ky? 
[-m 0? + (k4 = kp) — al = —k,Yi, 


(ka + ky — m,Q2)(kp — mgQ?) — kp? : 
a4 Eg = — ky Vi, 


kp an mp? 
A ao SS Se eo 
A (ka + kg — m,Q?)(kg — mgQ?) — kp? 


Hence 
2 eee ee 
(ka + kp — m4Q? (kg — mgQ?) — kp? 
and so, from Equation (S4) it follows that 

| = kek, Ysin Qt 


These results tally with the observation that both particles 
move with simple harmonic motion of angular frequency Q. 


XA 


Solutions to the exercises in Section 4 


; k 
1. Putting wo? = w,? = |— + — 
mz Mpg 


matrix equation (1), we have 


into the top line of the 


k k k k 
ee ta + X= 0, 


Ma my, Mg A 
XA XB 
so cx {), 
me Mg, 
Mpg 
1.e X4 = — > Xp 
Ma 


; ; cz m 
and so the required displacement ratio 1s — —. 
mM, 


2. The equations of motion are as follows. 


20(Xxz a X42) = 5X4 (particle A) 
—20(xp — X4) + 10(xc — Xp) = 6X%_ (particle B) 
——10(x¢ — Xz) = 4X. (particle C) 
3. The eigenvector equation is 
dd 4 0 
20 5 10 x x, 
6 6 Xp | = —w*| Xp 
10 10 X¢ X¢ 
4 4 
or alternatively, 
(-—4+ 7) 4 0 cS 0 
3.33 (—5 + w’) 1.67 Tj Le 0 
0 25 (—2.5+7)| | xc 0 


The characteristic equation of the model is obtained, as 
usual, by evaluating the determinant of the matrix on the 
left hand side and equating the result to zero. That is 


(—4 4+ w’)((—5 + w?)(—2.5 + a?) — (1.67 x 2.5)) 
— 4x 3.33 x (-2.54+ w7) =0 
so. (—4+ w7)(12.5 — 7.5m? + w* — 4.17) 
— 13.33 x (—2.5+ w?) =0 
or (w°® — 11.5w* + 38.33? — 33.33) 
— (—33.33 + 13.33m7) = 0 
ie. w° — 11.5w* + 25? =0. 


This is the characteristic equation of the model. 


4. First, we need to find the eigenveetors corresponding to 
each of the normal mode angular frequencies. We can do this 
by substituting each of the normal mode angular frequencies 
into the first and last of the three equations which constitute 
the matrix equation (2). 


First normal mode: w,* = 0.50. 
From the first equation 

— 6x4 + 4xz = —0.50x, 
ie. “4x, = 550x, 


a ee 
Ae 4 


From the third equation 
2.5Xp so 2.5X¢ = —0.50x- 
1.€. 2.5Xp = 2.00X¢ 


xe 32S 
2 XB 2.00 


Thus, putting x, = 1 we have Xp = 1.38 and 
Xc = 1.25 x 1.38 = 1.72, and so, the lowest normal mode 
eigenvector is [1 1.38 1.72]7. 
Second normal mode: w,? = 3.55. 
From the first equation 
— 6x, + 4x, = — 3.55x,4, 
iG Ax 2A5x 


2.45 
=o fo 
tS 


From the third equation 
Laks cee 25%e = — 3.55x¢ 
be 2S 105%. 


Xc 2 
— = —— = — 2,39. 
os, 
Thus, putting x, = 1, we have xz = 0.61 and 
Xc = —2.39 x 0.61 = — 1.47, and so the second normal mode 


eigenvector is [1 0.61 —1.47]7. 


Third normal mode: w3? = 9.46. 
From the first equation 
—6x, = 4Xxp = —9.45x4, 


1.e. 4xp = — 3.46x,, 
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From the third equation 
2.5Xp — 2.5xc = —9.45x¢, 
Le: 2.5Xp = — 6.96x¢, 


XC 25 
ef SE ee 
ee 


Thus, putting x, = 1, we have x, = —0.86 and | 
Xc = (—0.36) x (—0.86) = 0.31, and so the highest normal 
mode eigenvector is [1 —0.86 0.31]’. 


We can now write the general expression for the free motion 
of the three particles in terms of the eigenvectors and the - 
normal mode angular frequencies by substituting into 
Equation (4): 


oe 1 
Xp | =A, | 1.38 | sin(O.71t + ¢,) 
Xe L722 
l 
et 0.61} sin(1.88t + $3) 
— 1.47 
l 
+ A; |—0.86| sin(3.08t + ¢3). 
0.31 


Appendix 2: Solutions to the 
problems 


Solutions to the end of unit problems 


1. The given equations of motion are 


3 3 
= 74 — +(X4 — Xg) = mX4 


l 


fa 
=(i4 = Xp) = 7 = MXBgZ 


l 


where m = m, = mz = 0.05 and | = 0.2. 


For normal mode motion x, = —w?x, and X, = —w?x, so 
that 
T T 
2— — w? —_— 
ml ml Xa}. {0 
T T 2 XB = 0 : (SI) 
oe 2-——@ 
ml ml 
The characteristic equation for the model is therefore 
i @ 
2-——@? = 
ae 0 
ij ee ee 
-— 2— —@ 
ml ml 


T 
We can simplify this equation by putting p = ai to obtain 
m 


(2u — w*)(2u — w?) — p? =0, 
Le. 4? — 4um? + w* — p? =0, 
or w* — 4um? + 3y? = 0, 
sO. (w* — 3u)(w? — p) = 0. 
Hence, for the normal mode angular frequencies we have 


E ‘2 
WM; a 1.€. OO, = = 
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2 — 3 = : L. : = Ss 
WM yh Ro WM 
(a) Substituting m = 0.05 and / = 0.2 and putting w, = 300 
we obtain 


‘3 
0.05 x 0.2 


ie. “T= (300)? x 0.05 x 0.2 = 900, 


so the required tension is 900 newtons. 


= (300)’, 


(b) We know that w, = /30 so in this case 
@2 = ./3 x 300 = 519.62, 


so the larger normal mode angular frequency is about 
519.62rads7'. 


(c) To find the two normal mode displacement ratios we 
substitute each of the two values of normal mode angular 
frequency in one of the two linear equations which make up 
the matrix equation (S1). Using the top equation and 


; T 
starting with w,7 = a’ ¥ have 
m 


i & Tt T 
2— — —|xX4 — — Xz, = 0, 
ml ml ml 


T 
Le. —(X,4 — Xp) = 0. 
ml 
Hence xg = X4, So the first normal mode displacement ratio 
R, =1, or in eigenvector form [1 R,]’ = [1 1]’. 


This means that the particles will move in phase with equal 
amplitudes. 


. Pes ae 
Next, putting 2° = eat) into the same equation we obtain 
m 7 


Bora FE z 
2— — —]X4 — — Xp = 90 

ml mi ml 
cH 
le. ——(x,4 + Xg) = 0. 

ml 

Hence xg = —X,4, SO the second normal mode displacement 
ratio R, = —1, or in eigenvector form [1 R,]’ = [1 —1]'. 


So the particles again move with equal amplitudes but they 
will be phase-opposed. 


When the particles are furthest from their equilibrium 
positions their positions are as sketched below. 


Zz 
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We have just seen that 
w@, = 300 and 


with corresponding eigenvectors 


1 1 1 1 
Bight ak 
so from Equation (11) of Section 2 the general solution is 
x4 = C, cos 300t + D, sin 300t 
+ C,cos 519.62t + D, sin 519.62t 


xp = C, cos 300t + D, sin 300t 
— C,cos 519.62t — D, sin 519.62t. 


Ww, ~ 519.62 


(S2) 


We shall also need the velocities: 
X4 = 300(—C, sin 300t + D, cos 3001) 
+ 519.62(—C, sin 519.62t + D, cos 519.62t) 
Xp = 300(—C, sin 300t + D, cos 3001) 
+ 519.62(C, sin 519.62t — D, cos 519.62t). 


Substituting in the initial conditions x,(0) = 0, x,(0) = 1 and 
X4(0) = Xp(0) = O gives 


C,+C,=0 
C,-C,=1 
300D, + 519.62D, =0 
300D, — 519.62D, = 0 


which can be solved for C,, C,, D, and D,. From the 
second pair of equations we obtain 


D, = Dy = 0. 
From the first pair we obtain 
Pt 
2 2 


Substituting these values into (S2) gives the particular 
solution: 


x4 =4cos 300t — $cos 519.62 
Xp = 4cos 300t + $cos 519.62t. 
Thus when t = 2 the displacements of the particles are 
x4 =4cos 600 — $cos 1039.24 ~ —0.1 
xp = 4cos 600 + 4cos 1039.24 ~ —0.9 


or even more approximately x, = 0 (metres) and xz, = —1 
(metres). 


2. We know two of the eigenvalues: A = 0 (since this is a 
degenerate model) and 4 = —2. So, by inspection, we can 
write the characteristic equation in the form 


A(A +.2)(A2 + 44 + 2) =0. 


The other two eigenvalues are therefore obtained by solving 
the quadratic equation 4? + 44 + 2 = 0. We obtain 


—4+./16-—8 
A= a eee = —0.59 or —3.41. 
So the eigenvalues are: 
A, =0, A= —0.59, Az =-2, Ag= —3.41. 
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The normal mode angular frequencies are given by 


25k 
i= Sees That 1s 
m 
k k k 
wo; =0, wm, 207? j—, 0, = 1Ab fe eg le 
m m m 
(all in rads‘). 


To find the displacement ratios, we write down the equations 
of motion. The figure below shows the notation I shall use. 


The equations of motion are 


k(xp = Xa) am mx 4 
—k(xp — X4) + k(xc — Xp) = MXp 


k(xp — Xc) — K(Xc — Xp) = MX¢ 


—k(xp <i Xc) = MXp. 
Substituting X4 = —w?x,4, Xp = —w?xg, X¢ = —w*x,- and 
Xp = —W’Xp we obtain the normal mode eigenvector 
equations 
' mo)” 
Xp— Xa =e eA 
ze: 
ma)? 
(Xp —X4) + Xc = X= a XB 
ma 
itp — Xc) = he Hie k X¢ 
ma? 
{kp > Xe) a k Xp- 


Substituting the highest value of m into these equations we 
have: 


from the first equation 
Ap AS —3.41x,, 
$0: Xp = —ZAi ry, 


from the second equation 


—XBZ l 


t a1 


1 
| +Xc¢— Xp —3.41xp, 


1.€. = 1.41xp — Xp + Xe = —3.41xpz, 
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SO ~ Xo = —Xzg, 
= 241% ,: 
from the fourth equation 
—Xp + Xo = —3.41xp, 
i€. Xc = —ZA41 xp, 
SO Xp= —0.41x;- 


= (—0.41 x 2.41)x, 
= at, 


Hence putting x, = 1 we obtain the eigenvector 
[1 —2.41 2.41 —1]’ so the displacement ratios 
corresponding to the highest normal mode angular frequency 
Wa, are Rap = — 2.41, Rac = 2.41 and Rap = —{, 


3. (i) From the discussion on vibration absorbers in Section 
3 we know that the main mass will remain at rest when the 


k 
forcing angular frequency is —* where kz is the stiffness of 
m 


B 
the absorber springs and mz is the total absorber mass. 


5.9 x 10° = 
age 188rads~-°. 


When the speed of the shaft is such that the force has this 
angular frequency the main mass will remain stationary. 


Hence the angular frequency = 


(ii) When the main mass is stationary the total force on it 
must be zero. 


The gravitational forces on the main mass and on the 
absorber masses are cancelled out by the initial compression 
of the main springs; to put it another way, the gravitational 
forces are accounted for in determining the equilibrium 
position from which the vibration displacements are 
measured. It follows that when the machine is running, the 
force due to shaft unbalance must be cancelled out by the 
action of the vibration absorber. So the total maximum force 
which the absorber exerts on the main mass must be equal 
and opposite to the force that the unbalanced shaft exerts on 
the main mass. 


Now the amplitude of motion of each absorber mass is 
0.001 m — the maximum deflection of absorber spring from 
the equilibrium position. The stiffness of each absorber 
spring is 2.95 x 10°Nm_~’. It follows that the maximum 
force on the main mass due to each absorber spring is 


2.95 x 10° x 0.001 = 2.95 x 10°N. 


So the total maximum force on the main mass due to the 
absorber is 


5.9 x 10°N. 


And this must be equal in magnitude (though opposite in 
direction) to the force on the main mass due to the 
unbalance of the shaft. 
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